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1 Introduction 

During the last years the phenomenon of stabilization of atom in a strong laser field 
attracted a great deal of attention M . While the existence of the stabilization of atom has 
been clearly demonstrated in the numerical experiments the clear analytical criterion of 
stabilization is still absent. Usually it is assumed that stabilization condition is satisfied 
if the energy of the laser photon is larger than the electron coupling energy and the 
amplitude of electron oscillations in the field is large in comparison with Bohr radius § . 
However, the recent investigations of the corresponding classical problem demonstrated 
that stabilization remains also in the classical atom [Q, |4j , where the above conditions are 
violated. The physical explanation of this phenomenon and the condition of stabilization 
were given in ||] but the detailed explanation of the effect still remains an open problem. 
For a better understanding of this stabilization I introduce here a one-dimensional atom 
model which I will call Kramers model (having in mind that it arose from the Kramers 
- Henneberger transformation). Numerical analysis of this model allowed to construct 
an approximate Kramers map which describes the process of energy excitation and gives 
conditions of classical ionization. In some sense the obtained Kramers map is quite close 
to the Kepler map Q which describes the motion in the limit of relatively small field. 
Indeed, even in the strong field the change of the electron energy happens only when the 
electron passes near the nucleus while far from it the electron follows the Kepler orbit. 

The paper is constructed as follows. In the section 2 a brief description of the Kepler 
map is given since analogy with this map can be useful in the stabilization regime. In 
the section 3 a qualitative explanation of the stabilization is presented. The numerical 
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analysis of the introduced one-dimensional Kramers model and the derivation of the 
Kramers map are carried out in the section 4. In the section 5 I discuss analogy between 
the stabilization and the channeling of electrons in the crystal. In the conclusion the 
possibilities of experimental observation of stabilization of Rydberg atoms are discussed. 

2 Kepler Map 

After the pioneer experiments of Bayfield and Koch in 1974 || the problem of mi- 
crowave ionization of highly excited states of the hydrogen atom has been investigated by 
many groups (see [f?J and Refs. therein). The fast ionization observed in the experiments 
was really surprising since about 100 photons were required to ionize the atom. The typ- 
ical experimental conditions were the following: no ~ 70, eo = en^ ~ 0.05, ojq = ujuq ~ 1 
where no is the principal quantum number of initially excited state, e and u are the 
strength and the frequency of microwave field (here and below we use atomic units). The 
classical dynamics depends only on the rescaled values eo and ujq. 

For the understanding of the process of ionization in linearly polarized field it is 
convenient to use the one-dimensional atom model @, |§j, ||]. The investigations of one- 
dimensional model showed that for high microwave frequency (wn 3 > 1) the dynamics 
of the system, which originally is ruled by the continuous Hamiltonian equations, can be 
described by the Kepler map ||: 

N = N + ksm</>, 4> = (f) + 2iruj(-2ujNy 3/2 (1) 

Here k = 2.58e/u; 5/3 , N = E/uj has the meaning of the number of absorbed or emitted 
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photons {E is the energy of the electron), <p is the phase of microwave field at the moment 
when the electron passes near the nucleus. The bar denotes the new values of the variables 
after one orbital period. 

The physical reason due to which the motion can be quite accurately || described by 
the simple area-preserving map is the following: when the electron is far from the nucleus 
microwave field leads only to a small fast oscillations which doesn't modify the average 
energy and the Coulomb trajectory of the electron. The change of energy happens only at 
perihelion where the Coulomb singularity leads to a sharp increase of the electron velocity. 
Ionization takes place when the energy of the electron becomes positive after a pass near 
the nucleus N > 0. Then the electron goes to infinity and never returns back. Therefore 
for the map (1) ionization is equivalent to absorption of trajectories with N > 0. 

To find the chaos border in the Kepler map we can linearize the second equation in 
(1) near the resonant (integer) values of urn 3 obtaining the Chirikov standard map |J: 

N = N + ksm(f>, 4> = ( p + TN (2) 

with T = 67ro; 2 n 5 . The global chaos appears for K = kT > 1 that determines the critical 
field strength above which the classical atom is ionized. In this regime excitation goes 
in a diffusive way with the diffusion rate D = (AN) 2 / At = k 2 /2 where r measures the 
number of orbital periods of the electron. 

The first numerical and analytical investigations of the quantum one-dimensional atom 
model H showed that quantum effects leads to the suppression of classical diffusion. In- 
deed, in the quantum case the variables (iV, <f>) becomes the operators with the commuta- 
tion rule [N, <j)\ = —i and the system is locally described by the quantum kicked rotator 
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|1C]. The photon number is analogous to the level number in the kicked rotator and the 
excitation probability decreases exponentially with the number of absorbed photons so 
that the ionization rate is proportional to Wi ~ exp(—2Nj/l ( p). Here Nj = uq/2ujq is 
the number of photons required for ionization, = D = 3.33e 2 /a; 10//3 is the localization 
length. For l& « Nj quantum ionization is exponentially small in comparison with 
the classical value. However, for > iVj the diffusion is delocalized and the process of 
ionization is close to the classical one. 

In the 3-dimensional atom the Coulomb degeneracy leads to a slow motion along 
energy surface that allows to describe the excitation in energy also by the Kepler map with 
a small change of constant k. The motion along the energy surface has some additional 
integral of motion that explains the existence of localization in 3-dimensional atom || Q| . 



Recently the existence of localization in the 3d case was reconfirmed in [ 12 1 . 



Quantum localization of classical chaotic ionization has been observed in the mi- 



crowave experiments with hydrogen [|13|, 14] and rubidium [15] atoms. Numerical simula- 



tions with the quantum Kepler map [16] reproduce the 10% ionization threshold obtained 



in the laboratory ||13|. The theoretical prediction for the quantum derealization border 



was also observed in the skilful numerical simulations [12]. 

Being very successful in the description of energy excitation the Kepler map, however, 
cannot be applied for the case of very strong field. Indeed, in its derivation it was 
assumed that the change of energy after one kick kui is much larger than the energy 
of free oscillations e 2 /2oj 2 . This gives the condition of applicability of the Kepler map 
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picture ||: 

5a; 4 / 3 (3) 

Let us note that this condition is independent on the initial state since no doesn't enter 
directly in the expression for eati- 

In the one-dimensional case for e >> cati a collision with the nucleus, being unavoid- 
able, goes in a fast way like with an elastic wall leading to a prompt ionization H. In 
the two-dimensional case for zero magnetic quantum number m such collision also always 
takes place if the amplitude of free oscillations e/2uj 2 is larger than the unperturbed dis- 
tance between the electron and the nucleus in perihelion l 2 /2 (I is the orbital momentum). 
This gives the condition of prompt ionization for I > (3/oj) 1 / 3 [17]: 



e > uj 2 l 2 /A (4) 

where it was assumed that I is few times less than n. For I < (3/w) 1 / 3 ionization is ruled 
by the Kepler map and for eo > u>o 2 ^/2.6 prompt ionization takes place after one orbital 
period (see (1)). Therefore, there is no stabilization of classical atom in the strong field 
for m = 0. However, for high orbital momentum the atom remains stable up to very high 
field values. 



3 Stabilization Border 

While for the magnetic number m = ionization always takes place in a sufficiently 
strong field the case of nonzero m is much more interesting. Indeed, for the linear polar- 
ization of the field the projection m is an exact integral of motion and the created by it 
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centrifugal potential gives a possibility to avoid a collision with the nucleus. To analyze 
the motion in the strong field it is convenient to use the oscillating Kramers - Henneberger 
frame and cylindrical coordinates in which the Hamiltonian has the form: 

2 2 2p2 ( /9 2 + (z _^ sin(a;i) )2 ) i/ 2 

If the frequency of the nuclear oscillations is large enough (the condition will be given 
later) then in first approximation the nucleus can be considered as a charged thread with 
a linear charge density a slowly dependent on z: (t(z) = uJ 2 /(ire(l — (zuj 2 /e) 2 ) 1 / 2 ). Then, 
for small z and p the Hamiltonian of averaged motion takes the form ||, [|] : 

H «» = ? T + ? t + : $ + M*) In(^) (6) 

The constant under the logarithm takes into account that for p » e/co 2 the coupling 
energy becomes much less than uj 2 / e. From (6) one easily finds the position of the potential 
minimum p = ^'Kejlmjuj and the frequency of small oscillations Q = 2\f2uj 2 / (irem) (for 
z « e/uj 2 ). The depth of the potential or the energy required for ionization of atom is 
approximately I ~ 2uj 2 L/ire with L = ln(2e/ ' (eiru) 2 m 2 )) /2. The minimal distance between 
the nucleus and electron is determined by the condition / = m 2 /2(p m i n ) 2 giving: 

m Fne 

The physical reason for the growth of the minimal distance with the field strength is 
the following: with the increase of the field the amplitude of the field oscillations grows 
leading to the decrease of attractive Coulomb force while the centrifugal potential remains 
the same. 
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The averaged description of the motion (6) is correct if the frequency of field oscilla- 
tions ijj is much larger than the frequency Q of oscillations in p. In that case the averaged 
Hamiltonian (6) is the constant of the motion with adiabatic accuracy and ionization of 
atom doesn't take place. This gives the stabilization border || 

e > e stab = f3— (8) 
m 

where (3 is some numerical constant. The same estimate can be obtained from the con- 
dition that the change of energy AE during the collision between the electron and the 
nucleus is smaller than /. Indeed, the change of the momentum is Ap ~ At/p m i n 2 ~ 
uj/(ep m i n ) and the change of the energy AE u> 2 / (e 2 p m in 2 ) is less than / if (8) is 
satisfied. It is interesting to note that the stabilization border (8) can be written as 
v n = e/u > v max where v n is the typical velocity of the nucleus and v max = 2/m is the 
maximal velocity of the electron in the atom without the external field. 

Another condition intrinsically used in the derivation of (6) and (8) is p m in < e/^> 2 
which gives e > m 2 oj 2 . Also, there are two qualitatively different situations depending on 
the ratio between m and (3/u;) 1 / 3 . In the case m « (3/w) 1//3 (stabilized atom regime) 
we have m 2 u 2 « 5a; 4 / 3 << j3uj/m = e s t a b- For small field amplitude (3) the excitation 
is described by the Kepler map and the complete ionization after one orbital period of 
the electron takes place for eo > wo 2 / 3 /2.6 |J. Between this border and above chaos 
border e c o = l/49o;o 1 / 3 ionization goes in the diffusive way which is also relatively fast. 
However, for the more strong field (8), when the Kepler map picture is not valid (see 
(3)), atom becomes stable. The case of opposite inequality is less impressive. Indeed, 
for m >> (3/iij) 1 / 3 (stable atom regime) we have f3uj/m « 5w 4 / 3 << m 2 uj 2 and atom 
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remains stable (nonionized) up to e ~ m 2 uj 2 as it was in (4) (I ~ m). Above this value 
a significant portion (order of half) of atoms will remain stable since condition (8) is 
satisfied. Finally, ionization takes place only when the value of p m i n (7) becomes larger 
than the size of the atom 2no 2 and the electron cannot be captured in the stable region 
during the switching of the field. This gives the destabilization border 

l6Luj 2 n 2 . . 

tdestab ~ 5 \^) 

This border is also valid for the case m « (3/u;) 1 / 3 . Of course, in that case the stabiliza- 
tion can be observed only for the time of field switching T sw less or order of one orbital 
period of the electron. Otherwise a collision with nucleus will take place at field strength 
e < tstab and atom will be ionized. 

The results of numerical simulation of ionization process of system (5) are presented 
on the Fig.l. The stabilization probability W sta b = 1- — Wi on is given for different field 
strengths eo and frequencies ujq. The ionization probability Wi on was defined as the 
relative part of the trajectories with positive energies after field pulse. The initial distri- 
bution of 100 trajectories (25 for ojq = 1000.) corresponded to a quantum state with fixed 
spherical quantum numbers (fixed actions and equipartition in conjugated phases). The 
initial value of orbital momentum was equal to Z/rio=0.3 and its projection was equal to 
m/wo=0.25. The time of field switching (on/off) measured in the number of field periods 
was chosen to be equal to T sw = loq (one unperturbed orbital period of the electron). 
The pulse duration of the field was Ti nt = 500u;o (500 orbital periods). The data clearly 
demonstrate the stabilization of atom for field strength larger than some critical value. It 
is convenient to define the stabilization border as the field strength e sta (,o(20%) for which 



W s tab = 0.2. The dependence of e s t a bo(20%) on ujq, extracted from the data of Fig.l can 
be well fitted by the theoretical expression (8) with (3 = 12 in the wide frequency range 
(see Fig. 2 of @). This dependence continues up to loq = 1000 where we enter in the stable 
atom regime with m > (3/w) 1 / 3 and where stabilization disappears in agreement with 
above theoretical arguments (see Fig. 1). However, the stability of atom in that case is 
of the other nature than it was in [17] since the condition (4) is strongly violated. So, for 
such strong fields the stability of atom is based on the same physical grounds (8) as in 
the stabilized atom regime for m « {3/uj) 1 / 3 . The numerical check of the dependence 
of stabilization border e s t a 6o(20%) on m as well as the destabilization border (9) on ojq 
demonstrates good agreement with the theory (8)-(9) 0. 



4 Kramers Map 

It is important to stress that according to (8) the stabilization can take place even 
when the size of electron oscillations a = e/u 2 is much less than the unperturbed size of 
the atom Uq. An example of the motion in this case is presented on the Fig. 2. In such 
a case the electron follows the usual Kepler elliptic orbit and his energy (the size of the 
orbit) can be changed only during his fast passage near the nucleus. In this sense we can 
expect that the motion can be effectively described by some map analogous to the Kepler 
map. 

To construct such kind of map let's introduce a simplified one-dimensional Kramers 
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model given by the Hamiltonian 

^ = ^ + 4 3— (10) 

2 2p 2 ( p 2 + ^( sin7 + sin (^))^i/2 

This model is obtained from the Hamiltonian (5) by neglecting the changes of z and 

considering z = —e/u> 2 siwy as a constant. In other words electron collides with the line 

p = always at the same z value. The physical reason for that is the following. The 

collision of the electron with the nucleus is analogous to a collision of a fast heavy particle 

with the light electron. In such collision the change of energy (velocity) takes place mainly 

in the perpendicular p-direction, while the velocity in z-direction remains practically the 

same. According to this physical picture the model (10) mainly presents the changes in 

p-direction. In that sense it is quite different from the well known one-dimensional atom 

model of Eberly j!J which implicitly takes into account the change of energy (velocity) 

only in z-direction. Also in [18] the authors considered the velocity change only in z that 



has led them to a higher stabilization border than (8), while the estimate for p m in has 
been found correctly (see (7)). 

According to the analogy with the Kepler map we can expect that the change of the 
electron energy in the model (10) will take place only when the electron passes near the 
nuclear and that it will depend only on the phase of the field cp = uot at that moment. 
If also the size of the orbit is much lager than the size of the nucleus oscillations (a = 
e/u> 2 « no 2 ) then the change of the phase is given by the Kepler law and is the same as 
in (1). Basing on these arguments we can assume that the dynamics of energy excitation 
is governed by the Kramers map of the following form: 

E = E + Jh(4>), 4> = (f) + 2iru(-2E)- 3 / 2 (11) 
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with E = ujN, where as in (10) N is the photon number, the maximum change of the 
energy is given by a constant J and the unknown function of the kick h(<p) varies in the 
interval [-1,1]. 

To check the validity of this map I integrated the continues equations of motion for 
the model (10) and plotted the change of energy as a function of the field phase at the 
moment when the value of p took one of its minimal values (p p = 0). Such approach 
allows to find the kick function h((f>) the examples of which are presented on Figs. 3,4. 
The numerical results clearly demonstrate that the function h exists. However, it has 
a quite unusual property. Indeed, some values of 4> never appear (even if the number 
of periods was increased in 20 times). These values of (f> are approximately equal to 
7r + 7, 2ir — 7 and correspond to that values of the field at which the nucleus passes via 
the point of collision z = — e/oo 2 smj. A more close consideration of motion near these 
special 0-values shows that the electron remains during some small time interval (within 
corresponding phase interval Acj)) near the nucleus making one (Fig.3) or two (Fig.4) 
oscillations in p of very small amplitude so that the value of p remains practically (but 
not exactly) the same. This gives correspondently two (or three) values of the phase 
4> with the same change of AE since the value of E was determined in the aphelion. 
This of course puts the question about the derivation of the Kramers map in some other 
synonymous form. However, the main properties of the motion can be derived already 
from the approximate representation (11) where we will define the function h in the empty 
intervals by connecting the last points at the ends of the interval by straight line. 

Defined in such a way the Kramers map has the properties quite similar to the Kepler 
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map. Indeed, the function h is close to coscp and the approximate chaos border in (11) 
can be defined by the linearization of the second equation giving: 

K = 6irujJn 5 > 1 (12) 

where we used substitution E = — l/2n 2 . According to this criterion and in agreement 
with the numerical data the motion is chaotic for the cases of Figs. 3,4. If to introduce 
k = J/(jJ, which will give the number of absorbed photons after an orbital period, we will 
get the same formulas for the diffusion rate D = k 2 /2, the localization length (I = D) and 
the ionization time Tj on = Nf/D as in the Kepler map. In this sense the most important 
problem is the definition of the dependence of J on the parameters of the system. 

According to the results of the previous section the amplitude of the kick J must 
decrease exponentially with the increase of the stabilization (adiabatic) parameter S = 
em/to ~ u/Q (see (7), (8)). This expectation is in agreement with the results presented 
on Fig.5. Indeed, the exponential decrease of J with the field strength, and therefore 
stabilization, are evident. Let us at first discuss the properties of J for nonzero values of 
7. Even though the value of energy for the cases of Fig.5 was quite small nevertheless 
there is some dependence of J on energy. An example of such dependence is presented 
on the Fig. 6. We see that InJ depends on energy E approximately in a linear way and 
goes to a constant value for E = 0. This means that in the limit n 2 » e/uj 2 the value 
of J is independent from no- This result is consistent with the above arguments that the 
change of the energy takes place only in the small vicinity of the nucleus. However, in the 
difference from the Kepler map it is necessary to have quite strong inequality —aE << 1 
to neglect the dependence of J on E. We will try to explain this fact later. In this regime 
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of small energies the main change of the phase of the field between collisions (the second 
equation in (11)) is obviously given by the Kepler law. 

To determine the dependence of J(E = 0) on the parameters it is convenient to fix 
the stability parameter S that allows to eliminate the strong exponential dependence and 
to find the factor before the exponent. The numerical results are presented on Fig.7. The 
values of J(E = 0) were obtained from nonzero energies by linear extrapolation to E = 
(see Fig. 6). The numerical data clearly show that for the fixed S the value of J(E = 0) 
is independent on the frequency and is inversely proportional to m 2 . In principle the 
factor 1/m 2 gives simply the correct dimensionality however the independence on the 
other dimensionless parameter v = muj 1 / 3 is not so obvious. 

Combining all the obtained numerical results we can present the dependence of the 
kick amplitude J on the parameters for | E \ e/J 2 « 1 in the following form: 



where 51,2,3 are some functions weakly dependent on 7. For 7 = 0.6 we have from Figs. 
5-7 that gi ps 0.13, g 2 ~ 0.19, (73 ps 0.08. The numerical data for other values of 7 show 
that the fitting parameters vary not more than in 2 times for practically the whole interval 
of 7. For example, g x = 0.1 and 0.2, g 2 = 0.13 and 0.21, g 3 = 0.045 and 0.1 for 7 = 1.2 
and 0.3 correspondingly. 

To understand the numerically obtained formula (13) for J it is possible to make 
the following estimate. Taking the partial time derivative from the Hamiltonian (10) we 
obtain the expression for the change of the energy after one orbital period: 



J = 



gisin-y 



exp(-(# 2 - gzeE/uj 2 )em/u) 



(13) 




cos(?7 + 0)(sin7 + sin(?7 + 4>))drj 



p 3 / 2 (r]) 



(14) 
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where rj = u>t and in the denominator we neglected the term with e 2 /w 4 in comparison 
with p 2 . We can assume that near the nucleus the time dependence of p is the same as 
for a free electron with momentum m that gives p 2 (t) = pQ + v 2 t 2 where po is the minimal 
distance from the center and v is the velocity of the electron far from the center. For 
this free motion with the fixed momentum m we have the relation: pg = m 2 /v 2 . For 
the velocity it is possible to use the following expression: v 2 = uj 2 /Ce + 2E. Where the 
first term takes into account the fact that the energy must be measured in respect to the 
minimum of the effective potential (see (6)) and C is some unknown constant. It is easy 
to see that C determines the minimal distance p\ = Cem 2 /uj 2 for E = 0. In principal the 
value of C depends on 7. 

After substitution of all these expressions in (14) we obtain the following estimate 
02 

J 97 1m sm~fexp{-CS/(l + 2CeE/LU 2 )), h{4>) = cos{4>) , S = em/w(15) 

Of course, the presented derivation is not exact. However, it reproduce quite well the 
exponential dependence (13) (while the factor before the exponent is not in agreement 
with the dependence obtained from the numerical simulation). Indeed, numerically h{4>) 
has maxima near and tt. Comparison of (15) with (13) gives g2 = C and 53 = 2C 2 . The 
value of C can be defined directly from the numerical simulation of the one-dimensional 
Kramers model for different 7. The comparison of the (72 with C is presented on the Fig. 8 
showing the good agreement with the prediction. The ratios of the numerical value of 
gs (see above) to the theoretical value 2g\ are equal to 1.13, 1.1, 1.33 correspondingly 
for 7 = 0.3,0.6, 1.2 and are also in good agreement with the theoretical estimate. In the 
future estimates we will use the expression (13) with the theoretical substitution for 52 
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and g%. Let us also mention that for 7 = we get from (14) that 52 = 2C (from Fig. 5 
the ratio to the theoretical value is approximately 1.1) and h((p) = sin(2</>) that is quite 
close to the numerical data. Further theoretical analysis is required to obtain the factor 
before the exponent in (13). 

While still there are some unclear questions with the construction of the Kramers 
map the approximate consideration made above and the analogy with the Kepler map 
allow to understand the main properties of motion. If the number of photons required for 
the ionization is large then, as it was for the quantum Kepler map, it is possible to have 
diffusive excitation and quantum localization of chaos. However, due to high values of 
the frequency it is also quite easy to have a situation when one photon can already lead 
to ionization. In this case for k = J/uo < 1 the one-photon ionization rate (per unit of 
time) is given by the perturbation theory and as for the Kepler map (see ||) it is equal: 
J 2 



(16) 



For J > w approximately a half of probability is ionized after one orbital period ( in (11) 
as in the Kepler map the orbit is ionized if after a kick E > 0). From (16) it is clear 
that we may have long living states if the field is sufficiently strong. From the quantum 
view point one of the most interesting cases is the case of small m. In this case we need 
to make substitution m — > m + 1 since as it is well known the correct quasiclassical 
quantization leads to the appearance of the effective centrifugal potential even for zero 
orbital momentum. That gives the stabilization border e > 10a; for m = 0. 

Finally let us mention that in (11) we assumed that J is independent on energy. To 
take into account this dependence we need to put in the first equation J = J(E) and in 
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the second equation to add the phase shift A(p = dJ / dE f ((/>) with h{<p) = 
that way the map will remain canonical. 



-df(4>)/d(j). In 



5 Channeling Analogy 

Here I would like to discuss the analogy between the phenomenon of stabilization of 



atom in strong field and the channeling of particles in a crystal (see for example [19 
and Refs. there in). Let's consider the electron moving in the crystal with the velocity 
v ~ c = 137 (we will consider nonrelativistic case). Then in the frame of the moving 
electron its interaction with the protons in the crystal lattice will have approximately 
the form (5) if to take into account the interaction only with a nearest proton. On the 
grounds of that analogy we find that the effective distance between atoms in the crystal 
a and the velocity of the electron are equal to: 

a = -4, v = - ( 17 ) 

The frequency of perturbation is uj = v/a so that e = v 2 /a. Since in the crystal the 
distance between the atoms is approximately the same in all directions the analogy is valid 
for e/uj 2 > riQ 2 . The necessary condition of channeling is that the critical injection angle 6 
must be much less than one that implies: 6 w v±/v « l/(vy/a) « (w 4/3 /e) 3/2 << 1. This 
is the condition of unapplicability of the Kepler map (3). From the stabilization condition 
(8) it follows that channeling takes place for electrons with momentum m > 10/v. This 
is always satisfied for fast electrons with v ~ 137. The existence of channeling for very 
energetic electrons (that corresponds to strong field for stabilization problem) gives one 
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more evidence for existence of stabilization of atom in strong field in the regime when one 
photon frequency is larger then the ionization energy. 

6 Conclusion 

Basing on the Kramers map (11) and using its analogy with the Kepler map we 
obtained the estimate for the one-photon ionization rate (16). This ionization rate sharply 
decreases with the stabilization parameter S = em/uj. Such stabilization for excited states 
has some interesting advantages in comparison with the stabilization of atom in the ground 
state. Indeed, in this case stabilization can take place with e << 1 and lo << 1. This 
leads to a large energy difference SE between the excited states and the ground state. 
So, the energy in an exited state is approximately (e/u)) 2 /2 » 1 while the energy of 
the ground state remains as in the unperturbed atom (it is not the case for e, u » 1 
when the ground state is also stabilized since there the electron has the same energy of 
free oscillations) . Due to that in the case of Rydberg stabilization it is possible to have 
radiative transitions to the ground state with the radiation of X-ray photons. For the 
frequency of CO2 laser 00 « 1/300 (0.1 ev) and m=0 (or 1) the stabilization will take 
place for e « 1/30 (1.6 10 8 V/cm). The size of the atom will be larger than the size 
of the field oscillations a = e/co 2 for n > 40. According to (13) and (16) for the field e 
= 5 10 8 V/cm and n=60 the life time of the atom will be about 5 10 5 orbital periods 
or 10~ 6 seconds (we take for the estimate the case with 7 = 0.6). Of course, to obtain 
such states the time of field switching must be less than the time of orbital period as we 
discussed above. Since recently it was predicted that the Rydberg atoms can form long 
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living states (bands) in the solid state [^] (giving very high density of excited atoms) it 
will be interesting to consider a possibility of stabilization not only for a separate atom 
but also for such Rydberg solid state. 

I had started to be interested in the problem of microwave ionization of hydrogen 
atom in far cold 1980 during the first visit of Jeff Tennyson to the group of Boris Chirikov 
at Novosibirsk. Now, after many years of researches by different groups, this problem still 
continues to live by its own life as well as the memory about Jeff continues to live among 
the people who met him in Siberia. 
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Figure captions 

Fig. 1: Stabilization probability W s t a b = 1- — Wi on is given for different field strengths 
e and frequencies (w =0.3 (o), f. (*), 3. (+), 10. (O), 30. (A), 100. (O), 300. 
(A), 1000 (•)). 

Fig. 2: Example of trajectory for ujq = 300, eo = 20000 with initial I /no = 0.3 and 
m/riQ = 0.25; 10 5 field periods are shown. 

Fig. 3: Example of numerically obtained kick function h{4>) in Kramers map (11) for 
e = 3 10 4 , oj = 125, m=0.2, 7=0.6, £=-0.125 (so that effective n =2), J=l.l 10~ 4 . 
Near 200 orbital periods (points) are shown. 

Fig. 4: The same as Fig.3 with e = 410 4 , 7=1.2 and J = 5.810" 4 . 

Fig. 5: Dependence of the kick amplitude J in (11) on stabilization parameter S = em/uo 
for uj = 125, m = 0.2, 7 = (o); u = 1000, m = 0.1, 7 = (+); w = 125, m = 0.2, 
7 = 0.3 (open squares); u> = 125, m = 0.2, 7 = 0.6 (points); 00 = 125, m = 0.2, 
7 = 1.2 (full squares). For all cases £"=-0.125. Lines are drawn to adopt an eye. 

Fig. 6: Example of dependence of J on energy I E I for e = 2.2 10 4 , u = 125, m = 0.2, 
7 = 0.6 (points). 

Fig. 7: Dependence of J on m for fixed stabilization parameter S = 35.2 and £=0; 9 
cases are shown for uj in the interval [10,1000] and m in the interval [0.05,0.6]. The 
straight line shows the dependence J ~ 1/m 2 . 
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Fig. 8: Dependence of C = poe/S 2 on 7 (full line). Points give values of §2 to demon- 
strate theoretical relation C = <?2- 



23 



